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Abstract. We present new near-IR light curves for six Cepheids in the young blue LMC cluster NGC1866 as well as high 
precision radial velocity curves for ten Cepheids in NGC1866 and two in NGC2031. For the six Cepheids in NGC1866 with 
new J and K light curves we determine distances and absolute magnitudes by applying the near-IR surface brightness method. 
We find that the formal error estimates on the derived distances are underestimated by about a factor of two. We find excellent 
agreement between the absolute magnitudes for the low metallicity LMC Cepheids with the Period-Luminosity (P-L) relation 
determined by the near-IR surface brightness (ISB) method for Galactic Cepheids suggesting that the slope of the P-L relations 
for low metallicity and solar metallicity samples could be very similar in contrast to other recent findings. Still there appears 
to be significant disagreement between the observed slopes of the OGLE based apparent P-L relations in the LMC and the 
slopes derived from ISB analysis of Galactic Cepheids, and by inference for Magellanic Cloud Cepheids, indicating a possible 
intrinsic problem with the ISB method itself. Resolving this problem could reaffirm the P-L relation as the prime distance 
indicator applicable as well to metallicities significantly different from the LMC value. 



Key words. Cepheids; Magellanic Clouds; Stars: distances; 
Stars: fundamental parameters; Galaxies: distances and red- 
shifts 



1. Introduction 

Cepheids still provide one of the most important steps on the 
extragalactic distance ladder through the Period-Luminosity 
(P-L) relation. One distinguished example is the major effort to 
determine the value of the Hubble constant to better than 10% 
using the Hubble Space Telescope as presented by Freedman 
et al. ||2001 ). Still, the zero-point of the P-L relation remains 
an issue of heated debate and Freedman et al. (2001 ) chose to 
adopt a zero-point based on a best estimate of the distance to 
the Large Magellanic Cloud of 18.50+0.1 mag, and to assume a 
weak, but still ill-defined, effect on the zero-point as a function 



of metallicity. For a recent summary of the spread of current 
LMC distance estimates see e.g. Gibson J2000> and Benedict 
et al. (f2002t. 
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Fouque et al. (IMBl and Storm et al. ( I2004b> (S04 here- 
inafter) have determined direct distances and absolute magni- 
tudes to a large sample of Galactic Cepheids using the near-IR 
surface brightness (ISB) method as calibrated by Fouque and 
Gieren (I1997> . Unfortunately the slopes of the derived Galactic 
P-L relations differ from those observed in the LMC by OGLE 
(Udalski et al. .l999j in the optical and by Persson et al. (.2004J 
in the near-IR. If this effect is caused by the difference in metal- 
licity of the two Cepheid populations it will have harmful con- 
sequences for the use of the P-L relation as a distance indicator. 

Here we present new data for a sample of short period LMC 
cluster Cepheids and extend the ISB analysis to six Cepheids 
in the blue populous LMC cluster NGC 1866. This analysis will 
allow us to compare directly the absolute magnitude of LMC 
Cepheids with those for Galactic Cepheids using exactly the 
same procedure for both samples. It will also provide a first di- 
rect comparison of the slope of the P-L relations for low metal- 
licity stars in the Magellanic Clouds with the P-L relation for 
near solar metallicity Galactic Cepheids by including the SMC 
Cepheids studied by S04. 



2 



Jesper Storm et al.: The near-IR Surface Brightness Method applied to six Cepheids in the young LMC cluster NGC1866 



Additionally we will exploit the fact that the NGC1866 
Cepheids can all be considered to be at the same distance so we 
can measure directly the dispersion of the measured distances 
and thereby the true random error estimate for the method, 
something which has previously not been possible. 

NGC1866 is a unique object for many types of studies. It is 
young 100 Myr, Brocato et al. '1989'), and it is at the same 
time metal-poor when compared to similarly young Galactic 
stars. Hill et al. (|2000l find from high resolution spectroscopy 
of three giants a value of [Fe/H] = -0.5 + 0.1 confirming 
the photometric measurement of [Fe/H] - -0.46 + 0.18 by 
Hilker et al. J1995> . It is also very populous and thus provides 
an excellent laboratory for studying stellar evolution for inter- 
mediate mass stars (e.g. Brocato et al. ( 2003 1 and references 
therein). In the present context the most notable feature is the 
large population of cluster Cepheids. The possible association 
of Cepheids with NGC1866 was noticed early by Shapley and 
Nail (I1950> and photographic light curves were published by 
Arp and Thackeray ( 1967 1 for seven stars, six of which are an- 
alyzed in the present paper. Additional Cepheids were identi- 
fied by Storm et al. (19881, Welch et al. ([TWTl (W91 here- 
inafter), Welch and Stetson (1993) and most recently, using 
Hubble Space Telescope data, by Brocato et al. C2004 i. The 
Cepheids in the most crowded parts of the cluster are presently 
too crowded to be studied accurately from the ground, still the 
list of Cepheids presented by Welch and Stetson ( .1993 1 has 
more than 20 entries of which most have determined periods. 
NGC1866 thus provides us with a unique laboratory to study 
the Cepheid instability strip, which has already been exploited 
by e.g. W91, BertelU et al. (tT993t . and Brocato et al. ( l2U04l . 

Cote et al. (I1991> performed the first Baade-Wesselink type 
analysis of Cepheids in NGC1866 based on optical data, and 
the present analysis will build on their observational data as 
well as our new data. W91 provide an excellent historical sum- 
mary of the research on the interesting cluster NGC1866 up to 
the time of that paper 

Gieren et al. ( I2000bt presented a first near-IR surface- 
brightness analysis for the NGC1866 Cepheid HV12198, and 
we present here a re-analysis of this star based on the same pre- 
cepts as employed for the other Cepheids in the Milky Way and 
in the SMC. 



2. The Sample 

For the present work we have obtained J and K light curves 
as well as radial velocity curves for a sample of Cepheids in 
the young blue populous clusters NGC1866 and NGC2031. 
Finding charts for the NGC 1866 stars can be found in Storm 
et al. (119881 (for NGC1866-V4 see W91 though) and for the 
NGC 2031 stars in Mateo ( 1992 1. 

Only the data for the NGC 1866 Cepheids with Harvard 
Variable (HV) numbers proved of sufficient quality and com- 
pleteness to warrant the application of the infrared surface 
brightness method. 



3. The optical photometry 

Optical light curves using CCD detectors have been obtained 
for a number of NGC 1866 Cepheids by Walker ( 1987 1, W91, 
and most recently by Gieren et al. ( 2000a i. As shown by Gieren 
et al. (2000a I these data sets are all in very good agreement, 
especially in the V-band, and taken together they provide ex- 
cellent phase coverage and high accuracy for the six Cepheids 
which are of main interest to the present paper. Consequently 
we have adopted this combined BVRI data set for further anal- 
ysis. To ensure that all the data are on the same system we have 
applied small offsets (few hundredth of a mag) to the various 
data sets to bring them into optimal agreement with the Gieren 
et al. data set. 

The ephemerides which we have used are from Gieren et 
al. (2000a I for the stars with HV numbers, for the remaining 
stars we have determined periods and epochs on the basis of 
the available radial velocity data (see also Sec l5.3> and they 
are summarized in Tab|3] W91 found evidence for a slightly 
variable period for HV12198 but we do not find any significant 
variation over the 1 1 year period spanned by the optical data 
from Walker (1987i, W91, and Gieren et al. (UOOOai, so we 
have simply adopted a constant period for this star 

4. The near-IR data 

We have obtained near-IR photometry in the J and K' bands 
at ESO, La Silla, using the ESO-MPG 2.2m telescope with 
the IRAC-2 camera during observing runs in Dec. 1996 and 
Jan. 1998. We have also acquired photometry in the J and Ks 
bands at the Las Campanas Observatory (LCO) using the Im 
Swope telescope with the IRC in Nov. 1998 and Jan. 1999, and 
the 2.5m du Pont telescope in Jan. 1999. Finally we used the 
OSIRIS instrument at the 1.5m telescope at the Cerro Tololo 
Inter-American Observatory in Dec. 2000 to obtain / and Ks 
photometry. 

4.1. Data reduction 

For all the cameras we applied the appropriate non-linearity 
correction as quoted in the manuals. The exposure times were 
chosen to avoid exposing the detector to more than half of full 
well, so in all cases the applied non-linearity correction re- 
mained small. 

The data were obtained in sets where one set typically con- 
sisted of five or six different pointings on the target. In this way 
a local sky frame could be determined from the other frames 
in a set by median filtering. These local skyframes were sub- 
tracted from each of the original frames. In the case of more 
crowded fields, i.e. close to the cluster center, it was necessary 
to use sky frames based on data sets from another target with 
same exposure parameters obtained close in time. 

In case of the La Silla data the flat fields were determined 
by subtracting a domeflat with the quartz lamp off' from a dome 
flat with the quartz lamp on. A number of these exposure sets 
were averaged together to achieve a very high signal to noise 
ratio. A set of observations of a standard star placed in a grid of 
4 by 4 positions on the chip was used to determine the illumina- 
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Table 1. The K and { J — K) photometric measurements for the 
program stars transformed to the CIT system as described in 
the text. The complete table is available electronically from the 
CDS. 



ID 


HJD 


K 


0"(K) 


(J-K) 


0"(J-K) 




rdavsl 


[mas?! 


rmas?! 


fmas?! 


frnas^l 


HV12197 


2450438.65133 


14.524 


0.039 


0.481 


0.045 


HV12197 


2450439.57052 


14.428 


0.026 


0.327 


0.028 


HV12197 


2450439.65912 


14.387 


0.024 


0.354 


0.026 


HV12197 


2450439.76887 


14.399 


0.018 


0.362 


0.020 


HV12197 


2450439.86449 


14.408 


0.030 


0.365 


0.032 


HV12197 


2450440.55941 


14.382 


0.027 


0.417 


0.030 


HV12197 


2450440.67597 


14.387 


0.022 


0.427 


0.027 


HV12197 


2450440.79185 


14.402 


0.017 


0.458 


0.028 


HV12197 


2450441.76019 


14.610 


0.015 


0.403 


0.032 


cont. 













tion correction to be applied to the flatfield exposures and this 
corrected flatfield was used to flat field the science exposures. 

In the case of the Las Campanas data we used twilight flat- 
fields. Here the signal to noise ratio was somewhat lower, as the 
exposure level could not be as finely adjusted as was the case 
for the La Silla data. No additional correction for illumination 
effects was included. 

Finally the bad pixels were eliminated using the IRAF 
fixpix routine. 

4.2. Photometry 

The PSF-fitting photometry was carried out with DoPHOT2.0 
(Schechter et al.[T993} running on a DEC-Alpha with OSF/1 
UNIX. 

For each field a database was constructed containing the 
photometry of the stars from all the individual exposures of the 
field. On the basis of this database comparison stars could be 
selected and used as a reference to transform the photometry 
of each exposure to the instrument system. The detailed trans- 
formation procedure is described in the Appendix and the cal- 
ibrated photometric data is tabulated in Tab^and the resulting 
light and color curves have been plotted in Figs lll6l 

5. The radial velocities 

The radial velocity curves are based on high resolution spec- 
tra obtained with the Las Campanas Observatory 2.5m du 
Pont telescope equipped with the echelle spectrograph and 
the 2D-Frutti photon counting system. This system provides 
a dispersion of 0.2nm/mm (x 0.006nmpix"') and was used 
with a slit of 1 .5 x 4 arcsec which transforms into an instru- 
mental profile with a FWHM of approximately 3.5 pixels or 
12 km s ' corresponding to a resolving power of R ~ 20000. 

We also used the Cerro Tololo 4m Blanco telescope with 
the echelle spectrograph, the red long camera and the Tek2048 
CCD detector. This system provides a very similar resolving 
power {R w 20000) with a dispersion of about 0.25nm/mm 
which combined with 0.024 mm pixels again gives about 
0.006nmpix ^ The slit was again set at 1.5 arcsec width. 
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Fig. 1. K-band light curve and (V - K) and (J - K) colour 
curves for HV12197. 
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Fig. 2. K-band light curve and (V - K) and {J - K) colour 
curves for HV12198. 
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Fig. 3. K-band light curve and {V - K) and {J - K) colour 
curves for HV12199. 
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Fig. 4. 

curves 



K-band light curve and (V - K) and (7 - K) colour 
forHV12202. 
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Fig. 5. 

curves 



K-band light curve and (V - K) and (7 - K) colour 
for HV 12203. 




Fig. 6. K-band light curve and (V - K) and {J - K) colour 
curves for HV 12204. 



Integration times for our science objects with both systems 
were between 20 and 45 minutes depending on the weather 
conditions and the brightness of the star. 

5.1. The observations 

The Las Campanas data were obtained during three observing 
runs in Dec. 1995, Jan. 1996 and Jan. 2000, and the CTIO data 
were obtained during a single run in Jan. 2000. The correspond- 
ing Heliocentric Julian Dates are tabulated in Tab|2]with the 
derived velocities. 

The science spectra were interlaced with Th-Ar lamp spec- 
tra, and as the science objects were all in the same part of the 
sky, most science exposures had bracketing calibration spec- 
tra. A number of radial velocity standard stars were observed 
during each run to check the zero point. During the CTIO run, 
twilight sky exposures were obtained for the same purpose. The 
bright star AT67A (Arp and Thackerav. 1 1 9 67> . earlier selected 
by W91 to be similar to the Cepheids in spectral type and with 
sharp lines, which is located close to the two stellar clusters was 
observed frequently and served as our cross-correlation tem- 
plate star. 

5.2. Data reduction 

The Las Campanas data were first divided by a lamp flat 
field to take out the pixel to pixel variation in the detector. 
The apertures corresponding to each echelle order were ex- 
tracted using the IRAF' noao . imred . echelle . apall com- 
mand where the individual apertures had first been traced on a 
well exposed stack of standard star frames. The calibration Th- 
Ar apertures were extracted for each science exposure using 
the apertures previously defined by the corresponding science 
exposure. 

The wavelength solution was determined for each calibra- 
tion Th-Ar spectrum and using the IRAF refspectra and 
dispcor commands the science spectra were wavelength cali- 
brated by weighting the bracketing Th-Ar spectra according to 
their closeness in time to the science exposure. 

Ten echelle orders covering the wavelength range from 
444 - 5 1 3 nm were extracted and continuum corrected before 
they were merged into a one-dimensional spectrum using the 
scombine command. Finally outlying points, mostly from the 
low signal-to-noise ends of the spectral orders were removed to 
eliminate spurious signals. 

The CTIO CCD data were bias subtracted with the bias 
level corresponding to each individual quadrant of the quad- 
read CCD and each quadrant had its gain normalized. No pixel 
to pixel variation correction (flatfielding) was deemed neces- 
sary, and we preferred to refrain from such a correction to avoid 
the introduction of additional noise into the data. For the longer 
exposures it was necessary to remove the cosmic ray hits, and 
this was achieved by simply eliminating counts higher than a 



' IRAF is the Image Reduction and Analysis Facility, made 
available to the astronomical community by the National Optical 
Astronomy Observatories, which are operated by AURA, Inc., under 
cooperative agreement with the National Science Foundation. 
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Table 3. The ephemerides for the stars. 



ID 


Period 


Epoch of max. V light 




[days] 


[HJD] 


HV12197 


3.14381 


2450071.72 


HV12198 


3.52279 


2450069.59 


HV12199 


2.63918 


2450103.72 


HV 12202 


3.10112 


2450071.56 


HV 12203 


2.95411 


2450071.715 


HV 12204 


3.43876 


2450069.49 


NGC1866-V4 


3.31895 


2450068.5 


NGC1866-V6 


1.93386 


2450068.3 


NGC1866-V7 


3.388369 


2450067.9 


NGC1866-V8 


2.00714 


2450067.3 


NGC2031-V4 


3.42871 


2451572.7 


NGC2031-V11 


2.8244 


2451571.8 



certain threshold together with the pixels within a 2 pixel ra- 
dius. 

The spectral retrieval proceeded as for the Las Campanas 
data and 10 echelle orders covering the range 480 - 529 nm 
were finally extracted. The length of the extracted one- 
dimensional spectra were in both cases limited by the fact that 
the cross-correlation procedure used can only handle spectra 
with a maximum of 8192 pixels. 

The spectra were cross-correlated numerically using the 
technique developed by Tonry and Davis ( 1979) implemented 
in the rvsao package (Kurtz and Mink 1998, 1 within IRAF. 

For each run the spectra of the local reference star AT67A 
after dispersion correction were stacked to provide a high sig- 
nal to noise reference template spectrum. The science spectra 
were then cross-correlated against this template, for which we 
adopted a velocity of -0.7 km s"' as already found by W91 and 
confirmed by our measurements with respect to a twilight sky 
template and the standard stars observed. 

Nightly radial velocity offsets were determined on the basis 
of the shifts of the individual AT67A observations and/or twi- 
Ught sky spectra against the template and these offsets were fi- 
nally added to the derived velocities. Using the rvsao package 
the observed radial velocities were transformed into barycen- 
tric velocities, which are the values tabulated in Tab|3 The 
phases listed in the table are based on the ephemerides from 
TabE 

5.3. Radial velocity curves 

The radial velocity curves are shown in Figs l7l8l where our 
observations are shown with filled symbols and open symbols 
are used for data from the literature. We find excellent agree- 
ment with the data from W91 and with Storm et al. ( I2004a> 
for most of the stars. For HV12203 it appears that one of the 
observations from W91 is a misidentification (phase=0.68) and 
it has been rejected in the following. Two stars, HV 12202 and 
HV 12204 show a very systematic shift between observing runs, 
and we strongly suspect this is due to orbital motion. The case 
for HV 12202, as earlier noted by W91, is particularly strong 
as the offset is present in several observing runs. In the case 
of HV 12204 the shift is only seen for the W91 data. On the 



other hand we see no reason why the velocities for a single 
star should be systematically off when the rest of the obser- 
vations from that run are in good agreement with the other 
runs. To provide a radial velocity curve which is useful for a 
Baade-Wesselink analysis we have determined the shifts for 
HV12202 to be H-21 km s"' for the W91 data from the pe- 
riod 2447860 < HJD < 2447874 and our observations from 
the period 2451571 < HJD < 2451575 should be shifted by 
-H18 kms"' . In the case of HV 12204 the W91 data from the 
same period as for HV 12202 should be shifted by +7 kms"' . 
The resulting radial velocity curves are plotted in Fig|5] 
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Fig. 9. The radial velocity curves for the two NGC1866 
Cepheids with suspected orbital motion after the radial veloci- 
ties have been shifted as described in the text. The legend is as 
in Fig0 

We have determined the systemic velocity, Vy, for each 
star as the phase weighted mean velocity. These values are 
tabulated in Tab|3 For the two stars showing orbital motion, 
HV 12202 and HV 12204, we have performed the phase weight- 
ing on the original, unshifted, data. However, due to the limited 
sampling of the orbital motion a significant, but unknown, un- 
certainty remains on the systemic velocity for these stars. 

For the two Cepheids in NGC 2031 and for the four 
NGC1866 stars, V4, V6, V7, and V8, we have determined 
rough ephemerides as tabulated in Tab|3lon basis of the radial 
velocity curves. The epochs of maximum light have been esti- 
mated by forcing the minimum radial velocity to occur at phase 
0.0, which seems reasonable from the radial velocity curves for 
the other stars in Fig0and|9] 

6. The Infrared Surface brightness method 

The surface brightness method was developed in a series of 
papers by Barnes and Evans ( 1976 1 and Barnes et al. J 19761 
I1977t . They determined a calibration of the colour index (V-R) 
through a linear relation to the surface brightness parameter Fy 
defined as 

Fvi(p) = 4.2207 - Vo - 0.5 log 0((f)) (1) 
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Table 2. The barycentric radial velocity measurements, Vr, for each star, with the associated Heliocentric Julian Date, HJD, and 
phase, 0. The estimated random error for all the velocities is of the order 1km s ' . The complete table is available from the CDS. 



HJD 







HJD 


S 


* r 


HJD 





Vr 


HJD 


<l> 


Vr 


-2400000 






-2400000 






-2400000 






-2400000 






[days] 




[kms"' ] 






[kms"' ] 


[days] 


[kms"'] 


[days] 


[kms"' ] 


HV12197 






0.70 


320.1 


HV12204 




50096.5568 


0.58 


304.3 


D(JUOD.D04/ 


0.39 


298.5 


51 555 7082 


0.17 


286.3 


50066.5654 


0.15 


282.2 


50096.7672 


0.68 


306.5 




0.67 


314.2 


51 556 5402 


0.48 


308.8 


50066.7376 


0.20 


284.1 


50097.5561 


0.07 


293.2 


jUUo/.oU/I 


0.76 


315.5 


51556.6603 


0.53 


310.1 


50067.6675 


0.47 


302.3 


50097.7846 


0.19 


289.7 


50068.6029 


0.01 


283.8 


51556.8110 


0.58 


314.4 


50068.7943 


0.80 


318.9 








50095. 76js 


0.65 


310.9 


51557.6894 


0.92 


292.6 


50096.7144 


0.92 


304.9 


NGC2031-V4 




50096.6082 


0.92 


290.1 


51558.7151 


0.31 


297.4 


50097.7381 


0.21 


286.8 


51555.5 /9U 


U.Ul 


224.5 


50096.8013 


0.98 


281.6 


51572.7110 


0.61 


315.3 


51555.6398 


0.18 


284.3 


51555.7933 


0.07 


228.3 


50097.6318 


0.24 


292.6 


51574.7432 


0.38 


300.0 


51556.6832 


0.48 


304.0 


51556.5662 


0.29 


243.2 


50097.8436 


0.31 


292.1 








51558.5679 


0.03 


275.1 


51556.7349 


0.34 


245.8 


51555.6825 


0.03 


282.5 


HV12202 




51558.6887 


0.06 


276.0 


51557.7196 


0.63 


260.9 


51555.8249 


0.07 


278.6 


50066.6821 


0.43 


314.6 


51571.7364 


0.86 


323.3 


51558.5926 


0.89 
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0.10 


292.1 


50068.6640 


0.05 


293.5 


51572.7860 


0.03 


225.8 


51573.5842 


0.72 


316.6 


50097.6901 


0.43 


313.2 


50097.5985 


0.77 


312.6 


51573.5461 


0.25 


241.7 


51573.7776 


0.78 


315.2 


51571.6203 


0.72 


308.1 


50097.8082 


0.83 


306.7 


51573.6566 


0.28 


244.1 








51572.6177 


0.04 


275.1 








51574.5423 


0.54 


259.8 


HV12198 




51573.6939 


0.38 


293.0 


NGC1866-V6 




51574.6699 


0.57 


262.2 


50066.6203 


0.16 


284.0 


51574.6297 


0.69 


307.4 


50067.7598 


0.72 


306.2 


51574.7803 


0.61 


263.2 


50066.7786 


0.20 


289.4 








50068.6848 


0.20 


291.2 








50067.6057 


0.44 


304.1 


HV12203 




50095.7013 


0.15 


293.2 
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50068.7307 


0.76 
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0.27 
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50096.5415 


0.59 
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0.25 


236.7 
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0.44 


307.2 


50066.7585 


0.32 
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0.69 


311.4 
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0.32 


243.3 


50096.6464 


0.68 
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0.62 
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50097.5401 


0.10 
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51556.5884 


0.61 
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50096.8256 


0.73 


314.9 


50067.8480 


0.69 


319.6 


50097.7700 


0.22 


295.0 


51556.7571 


0.67 


265.4 


50097.6714 


0.97 


281.8 


50067.8637 


0.70 


319.2 








51558.6141 


0.33 


245.0 


51571.5604 


0.36 


301.0 


50068.7692 


0.00 


284.8 


NGC1866-V7 




51571.6610 


0.95 


221.2 








50095.8520 


0.17 


290.5 


50066.7189 


0.65 




51571.7738 


0.99 


Zi J.o 


HV12199 




50096.6908 


0.45 


308.3 


50067.6876 


0.94 


294.1 


51572.5861 


0.28 


241.3 


50066.6376 


0.95 


286.7 


50097.7193 


0.80 


318.3 


50068.6265 


0.21 


295.4 


51572.8168 


0.36 


244.6 


50066.8064 


0.01 


278.8 


51555.6186 


0.32 


300.9 


50096.5853 


0.47 


298.2 


51573.6204 


0.64 


263.1 


50067.5748 


0.30 


293.4 


51555.7429 


0.36 


303.3 


50096.7818 


0.52 


297.9 


51574.5781 


0.98 


220.5 


50067.7814 


0.38 


303.5 


51556.6311 


0.66 


319.4 


50097.5797 


0.76 


310.8 


51574.7004 


0.03 


222.1 


50068.5866 


0.69 


319.7 


51556.7839 


0.71 


322.0 


50097.8276 


0.83 


312.2 


51574.8164 


0.07 


223.5 


50068.8237 


0.78 


316.9 


51557.6672 


0.01 
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0.31 


295.9 
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0.39 


299.4 








50068.6428 


0.67 


311.7 









where is the stellar angular diameter, Vq the de-reddened vi- 
sual magnitude, and (p is the phase. 

The angular diameter is geometrically related to the dis- 
tance and stellar radius through 

e{(l>) = 2R(<p)/d = 2(Ro + AR((P))/d (2) 

where d is the distance and R is the radius. The variation of 
the radius can be determined by integrating the radial velocity 
curve of the star: 

A/?(0)= r -p[VA<p)-Vy]d(P (3) 



where p is the so called projection factor converting radial ve- 
locity into pulsational velocity, Vr{4>) is the radial velocity and 
Vy is the systemic velocity. 

As Fy is calibrated from interferometric measurements of 
stellar angular diameters it is then possible to solve Eq|2]for the 
two remaining unknowns, namely the distance, d, and the stel- 
lar mean radius, Rq. As discussed by Barnes et al. (2T)03) and 
S04, it is sometimes necessary to apply a small phase shift, A<p, 
to the radial velocity data to properly phase the spectroscopic 
and photometric data before solving Eq|2| 

Welch (11994') demonstrated the advantages of using the 
near-infrared colour index (V - K) together with the K-hand 
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Fig. 7. The radial velocity curves for six NGC1866 Cepheids. Our new data are marked with filled symbols, while data from the 
literature have been over-plotted as open symbols. The data are marked in the following way: LCO 1995: filled squares; LCO 
1996: filled diamonds; LCO 2000: filled circles; CTIO: filled triangles; Storm et al. ( 2004a.) : open circles; Welch et al. a991i) 



HJD < 2447520: open squares; Welch et al. ( [19911 2447860 < HJD < 2447877: open diamonds. 



magnitudes and made a first calibration of versus (V - K) 
on the basis of interferometric diameters of non-pulsating stars. 



Fouque and Gieren ( I1997t reformulated the near-infrared re- 
lation as a F\f, (y - K) relation which they calibrated, using 
an expanded sample of accurate interferometric measurements 



of non-pulsating giant stars. Later Nordgren et al. fl2002t and 
Kervella et al. C2004a) obtained very accurate interferomet- 
ric radii of Cepheids allowing them to calibrate the surface 
brightness method directly. For the Cepheid { Car Kervella et 
al. f 2004b) made a phase by phase comparison between their 
interferometric measurements and the values returned by the 



Fouque and Gieren relation. Recently Groenewegen ( I2004> has 
expanded and improved on the work on the non-pulsating stars. 
As there is excellent agreement among all these various cali- 
brations, we choose for consistency with S04 to still adopt the 
relationship from Fouque and Gieren ("1997'). In this way the 
left hand side of Eq|2]seems well established. 

The right hand side of the equation contains the radius vari- 
ation which, apart from the direct observables, depends through 



Eq|3lon one non-observable, namely the projection factor p. 
As an error in p carries directly over in the derived distance, 
it is of fundamental importance that p is well understood. S04 
and references therein argue for their choice of a weakly pe- 
riod dependent p, namely p - 1.39 - 0.03 logP where P is the 
pulsation period in days. 

6.1. Results 

Before applying the ISB method to the observational data we 
need to adopt reddenings for the individual stars. The distances 
derived from the ISB method are themselves very robust to er- 
rors in the assumed reddening, but the absolute magnitudes of 
course reflect directly any errors so introduced. 

We have adopted the canonical value of E(B - V) - 0.06 
which was originally determined by van den Bergh and Hagen 



d 19681 1 from integrated UBV photometry, and by Walker ( I1974> 
from stellar UBV photometry. We assume an uncertainty of 
about 0.01 mag on this result. 
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Fig. 8. The radial velocity curves for four additional NGC1866 Cepheids and for two NGC2031 Cepheids. The legend is as in 
FigEl 



Table 4. The systemic velocity, Vy, computed as the phase 
weighted mean velocity. The random uncertainty is of the order 
0.2 km s"' and the systematic error of the order 1 km s"' . 



ID 






[km s"' ] 


HV12197 


298.6 


HV12198 


299.7 


HV12199 


300.4 


HV 12202 


303.1" 


HV 12203 


302.4 


HV 12204 


296.5" 


NGC1866-V4 


302.6 


NGC1866-V6 


298.9 


NGC1866-V7 


298.7 


NGC1866-V8 


298.2 


NGC2031-V4 


247.7 


NGC2031-V11 


242.1 


a: Stars with suspected orbital motion. 



In Fig[TO|we have plotted the P-L relations in V, I, W, and 
K using the adopted reddening of £{8 - V) = 0.06. Over- 



plotted the data are the observed P-L relations based on the 
OGLE data (Udalski et al. I1999> as determined by Fouque et 
al. ( I2003> based on an average LMC reddening of E(B -V) = 
0.10 mag. The estimated relative random errors on the mean 
magnitudes in the optical bands are 0.005 mag and thus smaller 
than the points in the diagram. In the K band where the zero- 
points are determined individually for each star the estimated 
random errors are somewhat larger, of the order 0.04 mag. 

Groenewegen and Salaris (I2003> have recently claimed a 
reddening of E(B - V) = 0.12 towards NGC1866, but we do 
not see any support for this value in our own observational data. 
From Fig^Jit is clear that the V and / data appear slightly too 
faint (disregarding the slightly peculiar star HV 12204) but an 
increase of the reddening of only 0.01 mag would make the data 
in all the diagrams agree very well with the OGLE based rela- 
tions. Adopting an NGC1866 reddening of E{B - V) = 0.12 on 
the other hand, would cause all the data points in the reddening 
sensitive {V and /) diagrams to appear strongly (0.15 mag in V) 
over-luminous with respect to the OGLE based relations. We 
thus argue that the NGC1866 reddening of E(B - V) = 0.06 is 
in good agreement with the reddening scale adopted by Fouque 
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et al. (I2003I I which will also be used in the following. We stress 
the importance of using self-consistent reddenings and the su- 
periority of reddening insensitive indices hke W and better still 
the K band. 

On the basis of the observational data presented in the pre- 
vious sections we have applied the infrared surface brightness 
method, using the Fouque and Gieren J 19971 1 calibration fol- 
lowing the procedure explained in S04. The derived distance 
moduli and associated formal errors for the individual stars 
have been tabulated in Tab|5l together with the derived abso- 
lute magnitudes. We have also tabulated the phase shifts, A(p, 
which have been applied to the radial velocity data to match the 
photometric angular diameter curve as discussed by S04. 

For HV12198 we find a distance modulus of 18.31, which 
is shorter by 0.11 mag compared to the value of 18.42 + 0.1 
found by Gieren et al. ( 2000b i. This is largely due to the fact 
that, as discussed in detail by S04, to avoid effects from possi- 
ble shocks in the stellar atmosphere close to minimum radius, 
we only consider the phase interval between 0.0 and 0.8 here. 
An additional, but less important, reason is that we now use the 
bi-sector for fitting Eq|2l 

In Fig [H] we show as an example the fit to the data for 
HV12198 as well as the angular diameter variation as a func- 
tion of phase. 




Phase 

Fig. 11. Panel a shows the observed data set of angular di- 
ameters, 6 (in micro arc-seconds), versus radius variation, AR, 
corresponding to Eq|5|for the star HV12198. The straight line 
is the bi-sector fit to the data in the phase interval [0.0,0.8]. 
Panel b shows the angular diameter as a function of phase. The 
points correspond to the photometrically determined angular 
diameters and the curve coiTesponds to the integrated pulsa- 
tional velocity curve. In both panels the crosses represent dis- 
carded points. 
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Fig. 10. Panel a shows the observed intensity averaged K mag- 
nitude vs logP with the Persson et al. (20041 relation over- 
plotted. Panels btod show the observed V,I, and W magnitudes 
with the coiTesponding OGLE based relations from Fouque et 
al. (l2003t over-plotted. 
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Fig. 12. Panels a to c show the absolute magnitudes of the 
Cepheids versus log P with the Galactic relation from S04 over- 
plotted. Panel d shows the derived distance moduli as a func- 
tion of period. 



In Fig^lwe have plotted the derived absolute magnitudes 
in V and K as well as the Wesenheit index and the distance 
modulus as a function of log P with the Galactic relations over- 
plotted. It can be seen that the absolute magnitudes are in good 
agreement with the Galactic Cepheids, confirming the finding 
by S04 from SMC Cepheids that a metallicity effect on the ab- 
solute magnitude of these stars is small, if at all present. The 
Galactic relations are largely based on the data presented in 
S04 but expanded by four additional stars and some additional 
observational data from the literature, see Gieren et al. (i2005< 



for details. We note that the scatter in Fig^jis largely domi- 
nated by errors in the distance estimates as these correlate well 
with the deviations seen in the luminosities. The effect of color 
(temperature) on the luminosity due to the finite width of the 
instability strip is comparatively minor, as can be seen from 
FigEI which is distance independent and shows a very low 
scatter with only HV12204 exhibiting a significant deviation. 
We will later argue that this deviation can be explained by a 
blue companion. 
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Table 5. The measured distance moduli, radii, and absolute intensity averaged magnitudes for the NGC1866 Cepheids. The 
estimated error is just the formal error on the fit. The adopted reddening, E{B - V), and the phase shift, A(p, applied to the radial 
velocity data are also given. 



ID 


\ogP 


(m - M)o 




R 


0"fl 


Mb 


Mv 


M, 


Mj 


Mk 


Mw 


E(B - V) 


&(fi 






[mag] 


[mag] 


[«o] 


[Re] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 




HV12199 


0.421469 


18.336 


0.094 


25.0 


1.1 


-1.737 


-2.269 


-2.870 


-3.075 


-3.672 


-3.777 


0.060 


0.025 


HV 12203 


0.470427 


18.481 


0.092 


28.3 


1.2 


-1.989 


-2.552 


-3.153 


-3.562 


-3.930 


-4.060 


0.060 


0.050 


HV 12202 


0.491519 


18.289 


0.072 


28.5 


1.0 


-1.815 


-2.425 


-3.050 


-3.553 


-3.922 


-3.992 


0.060 


0.025 


HV12197 


0.497456 


18.165 


0.058 


25.9 


0.7 


-1.712 


-2.273 


-2.918 


-3.338 


-3.728 


-3.891 


0.060 


-0.020 


HV 12204 


0.536402 


18.202 


0.044 


28.3 


0.6 


-2.243 


-2.702 


-3.239 


-3.626 


-3.981 


-4.050 


0.060 


0.010 


HV12198 


0.546887 


18.314 


0.028 


29.8 


0.4 


-1.989 


-2.565 


-3.202 


-3.675 


-4.030 


-4.165 


0.060 


0.015 



7. Discussion 

7.1. Cluster membership of HV12204 
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Fig. 13. The colour-magnitude diagram for the NGC1866 
Cepheids in V, (B - V) and K, (V - K). The bluest object is 
HV 12204. 

Fischer et al. ( I1992> measured the velocity dispersion of 
NGC1866 from non- variable stars to be 2.3 kms ' and found 
a systemic radial velocity for the cluster of yy,NGCi866 = 
301 kms"'. W91 found a systemic velocity of Vy = 
293 kms ' for HV12204 and consequently they suspected it 
to be a non-member due to the inconsistent radial velocity and 
later works have all dismissed the star on this basis. However, 
with the additional radial velocity data we have presented here, 
it appears that HV 12204 is a single-lined spectroscopic binary 
exhibiting orbital motion. We do not have enough data to at- 
tempt to determine an orbit, but our data suggest that HV 12204 
could very well be a member of NGC1866. In Fig[TO]the de- 
reddened magnitudes of the stars have been plotted against 
logf with the OGLE based P-L relations from Fouque et al. 
J200 3^ and from Persson et al. f2004) over-plotted. It is clear 
from these plots that the NGC1866 Cepheids closely follow 
these relations. In the V and / bands HV12204 appears slightly 
over-luminous, but the effect does not seem to be caused by a 
difference in distance with respect to the other NGC1866 stars 
as the effect changes with wavelength and has disappeared in 
the K band as well as in the Wesenheit index. Consequently we 
do not find sufficient evidence to discard HV 12204 as a possi- 



ble non-member of the cluster. In Fig[21the colour-magnitude 
diagram for the NGC1866 Cepheids has been plotted and it can 
be seen that HV12204 is the bluest and brightest of the stars. 
We cannot exclude that the photometry is affected by an unre- 
solved blue companion which is further supported by the fact 
that the (J - K) color is in good agreement with that observed 
for the other cluster Cepheids. 

7.2. The distance to NGC1866 and the estimated 
distance errors 

We now have a sample of six NGC1866 Cepheids which we 
can use in the following to determine the distance to NGC1866 
and to the LMC, and to determine the size of the internal ran- 
dom errors of the method. 

The six Cepheids are all within a radius of 3 arcmin from 
the cluster center and assuming conservatively a maximum 
cluster radius of 10 arcmin, this corresponds to a physical ra- 
dius of about 150 pc at a distance of 50 kpc. Thus the maximum 
error due to depth effects within NGC1866 would be +150 pc 
on a total distance of 50 kpc which is +0.3 % or ±0.007 mag. 
For our purposes we can thus safely assume the stars to be at 
the same distance. 

We find a weighted mean distance modulus of 18.28 ± 0.05 
with a standard deviation of 0.11 mag. When we compare the 
standard deviation with the formal errors in Tab|5] which range 
from 0.028 to 0.094, it is obvious that the formal errors sig- 
nificantly (by about a factor of two) underestimate the true 
errors. It can also be seen from the table that the formal er- 
rors spread significantly, even though the data quality is rea- 
sonably uniform. In fact the variation of the estimated errors 
mostly reflect the overall agreement between the shape of the 
two angular diameter curves and to a much lesser extent the 
random noise of the data. In this way a few stars carry a pos- 
sibly unrealistic high relative weight. We thus prefer to use the 
unweighted mean as the most unbiased estimator, leading to 
(m - M)o = 1 8 .30 + 0.05^, in apparent excellent agreement with 
the main-sequence fitting distance from Walker et al. 120011 . 

Barnes et al. ("2003) have performed a very thorough statis- 
tical analysis of the optical surface brightness method using 
Bayesian techniques. A similar investigation of the infrared 
surface brightness method is currently in press (Barnes et al. 
.2005 ). Preliminary results indeed suggest that the formal er- 



- This is our best distance estimate at present, but we caution the 
reader to take into account the caveats in subsection l7.4l 
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rors of the type reported here are underestimates but that the 
derived individual distances remain unchanged. 



7.3. The LMC distance 

Even though NGC1866 is some way away from the LMC cen- 
ter we can attempt, at least formally, to determine the distance 
to the LMC on the basis of these data. 

As NGC1866 is a young object it must be close to the 
LMC disk. In Fig[TO|it seems as if the NGC1866 Cepheids are 
slightly (of the order 0.05 mag) fainter than the average LMC 
Cepheids in the V and / bands when disregarding HV 12204, 
but very close to the mean value in the reddening insensitive 
W index and in the K band. As previously discussed the de- 
viation in the V and / bands would disappear if the reddening 
were increased by AE{B - V) = 0.01, well within the esti- 
mated uncertainty of the reddening. Groenewegen and Salaris 



J2003> have suggested that NGC1866 is more distant than the 
surrounding disk population. At the same time the disk model 
of van der Marel and Cioni ( 1200 It at this position in the north 
eastern part of the LMC (RA: 05:13:39, DEC: -65:27.9 
(2000 . 0)) says that the disk is closer than the LMC center by 
0.059 mag. These two effects largely cancel out and confirm 
our finding that the NGC 1 866 Cepheids are close to the average 
LMC Cepheid distance, so we choose for the present purposes 
not to apply any correction for geometry, but add in quadrature 
an additional uncertainty of 0.05 mag. We thus formally find an 
LMC distance modulus of (m - M)o,lmc = 18.30 + 0.07 mag. 
We note that the derived modulus is on the short side of the 
available estimates, but as the following discussion will show, 
this should only be considered a preliminary estimate of the 
LMC distance modulus. 
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Fig. 14. The P-L relation for Galactic Cepheids with the 
NGC 1866 Cepheids overplotted. The dashed line is the OGLE 
P-L relation for an adopted LMC distance of (m - M)o = 
18.30 mag. 
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Fig. 15. The P-L relation for Galactic Cepheids with the 
NGC 1866 Cepheids overplotted. The dashed line is the Persson 
et al. (2004) P-L relation for an adopted LMC distance of 
(ot-M)o = 18.30 mag. 
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Fig. 16. The P-L relation for Galactic Cepheids with the 
NGC 1866 Cepheids overplotted. The dashed line is the OGLE 
P-L relation for an adopted LMC distance of (m - M)o - 
18.30 mag. 



7.4. Comparing the Period-Luminosity relations 

The absolute magnitudes which we have derived for the LMC 
Cepheids appear in excellent agreement with the values ex- 
pected from the Galactic relations derived using exactly the 
same precepts by S04. In figures I14I16I we have overplotted 
these stars on the Galactic Cepheids. We have also overplot- 
ted the SMC Cepheids from S04, which have a metallicity of 
[Fe/H] - -0.7 which is not very different from the metallic- 
ity of the NGC 1866 Cepheids of [Fe/H] = -0.5 from Hill et 
al. f2000''(. The filled lines in these plots represents the linear 
Galactic relations fit to the Galactic Cepheids. It is clear that the 
low metallicity stars follow this relation quite well and there is 
no obvious sign of a difference in slope between the low metal- 
licity and high metallicity populations. From the perspective of 
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the ISB method the two populations appear to follow similar 
P-L relations. 

However, the dashed lines overplotted in Figs ll4ll'Sl repre- 
sent the observed apparent P-L relations shifted to the formal 
LMC distance of (m - M)q - 18.30 mag. As already noted and 
discussed by Fouque et al. (2003 1, Tammann et al. (2003 1, and 
S04 these slopes differ from the slopes found from the Galactic 
sample as determined from the ISB method. Now we can see 
indications that these slopes also might differ from the slopes 
determined directly for the low metallicity stars from the ISB 
method. This suggests that there could be a problem either with 
the slopes of the apparent P-L relations, which is highly un- 
likely, or with the slopes derived from the ISB method. To clar- 
ify this point it will be necessary to study a significantly larger 
sample of LMC Cepheids spanning a large range of periods. 
Gieren et al. ( .2005J will present the first results from such an 
expanded sample. 

If we can indeed confirm the suggestion that the LMC and 
SMC Cepheids follow the same P-L relation from the ISB 
method as the Cepheids in the Galaxy, this will mean that the 
P-L relation is afterall a good distance indicator as the slope is 
not dependent on metallicity. However, the ISB method itself 
will then need mending and the distances and absolute magni- 
tudes derived here and in S04 will need revision. 

It is relevant to note that the theoretical models of Caputo 
et al. ( I2000t predict that the metallicity effect on the luminos- 
ity for short period Cepheids is very small, exactly as seen in 
Fig^S The predictions for the longer period stars, and thus for 
the slope, are tested with the expanded sample of Gieren et al. 
ll2005l . 

8. Conclusions 

We have presented new near-IR light curves for six NGC1866 
Cepheids and high precision radial velocity curves for six 
NGC1866 Cepheids and two Cepheids in NGC2031. 

Using these data we have shown that HV 12204 is most 
likely a member of NGC1866 but subject to orbital motion, 
just likeHV12202. 

Following the precepts in S04 we have derived absolute 
magnitudes and distances to the six NGC1866 Cepheids. The 
distances exhibit a standard deviation of 0. 11 mag suggesting 
that the formal eiTor estimates on the individual distance mod- 
uli are underestimated by about a factor of two. 

Formally we find a distance modulus of (m-M)o = 18.30 + 
0.05 to NGC1866, but due to the potential problems with the 
method in its present form uncovered here, this should be re- 
garded as a preliminary estimate. Similarly the formal LMC 
modulus of (m - M)o = 18.30 + 0.07 might need revision. 

The absolute magnitudes My, Mk and the Wesenheit in- 
dex Mw for the NGC1866 Cepheids are all in excellent agree- 
ment with the P-L relations determined by S04 for Galactic 
Cepheids. 

The NGC1866 Cepheids with a metallicity of [Fe/H] = 
-0.5 and the SMC Cepheids from S04 with a metallicity of 
[Fe/H] = -0.7 appear to follow largely the same P-L relation 
as the Galactic Cepheids, in particular they seem to follow a 
similar slope. However, this slope is different from the slope 



of the directly observed apparent P-L relations in K, V, and W 
suggesting an intrinsic systematic eiTor with the ISB method 
which casts doubt on the distances and absolute magnitudes 
derived by the method. On the other hand if this problem orig- 
inates with the ISB method itself, this might resolve the very 
unpleasant problem of the slope of the P-L relation depending 
on metallicity. 
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Appendix A: Photometric transformations 

A.1. Transformation to the common instrument system 

The JK data from all the observing runs were transformed to 
the Las Campanas 1998 instrument system before being trans- 
formed to the standard system. 

A.1 .1 . The Las Campanas and CTIO data 

In the case of the Las Campanas observing runs and the 
CTIO run, the colour terms for the transformations to the Las 
Campanas instrument system were determined using a data set 
for the LMC cluster NGC 2136. This field contains a large 
number of stars without becoming seriously crowded. The cen- 
tral part containing the cluster itself was disregarded to reduce 
the number of stars with systematic errors due to poor sky sub- 
traction or neighbour subtraction. For each set of J and K ex- 
posures the transformation to the Las Campanas (LCO) instru- 
ment system was determined, assuming relations of the form 

^^LCO inst = ^^inst + aK( J - K)in,t + JSk (A. 1) 

^LCO inst - Jimt + aj( J - K)mst + j3j (A.2) 

For each observing run the resulting slopes (aj and a^) 
were averaged. The Las Campanas data from the other observ- 
ing runs did, as expected, not exhibit significant colour terms, 
and we have assumed a value of zero for a in the actual trans- 
formations. For the CTIO data there were slight colour terms 
of UK = -0.03 + 0.01 and aj = -h0.03 + 0.01. 

The individual J and K data sets were then reanalyzed with 
the slopes fixed at the adopted values leaving only the zero 
point offsets p j and Pk to be fitted. Only stars with estimated 
eiTors less than 0.05 mag were used in this fit and the known 
variables were of course excluded. 
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A.1.2. The ESOdata 

NGC2136 was not observed during the ESO observing runs. 
Instead we observed 1 1 standard stars from the list of Carter 
and Meadows (1995), during the night Dec.21, 1996, to de- 
termine the colour terms. The stars span a colour range from 
(J - K) = -0.02 to 0.96. For the /T-band we found no signifi- 
cant colour term and in the 7-band we found ay - -1-0.06+0.02. 
This is in good agreement with the value quoted in the user 
manual for this colour range. As the colour range for our sci- 
ence targets as well as for the comparison stars in these fields 
are all very similar (±0.2 mag), the exact value of the colour 
term is not very important, as long as it is small. 

Carter (1993) gives the transformation for normal stars to 
the CIT/CTIO system as 

Ken ^ Kc- 0.021(7 - K)c - 0.002 (A.3) 
(J - K)cn = (7 - K)c - 0.105(7 - K)c - 0.001 (A.4) 

with standard errors of 0.01 mag. Inserting the previously 
found colour terms for the ESO data in these equations gives: 



(J - K)c 


= (7 -/:)eso + 0.06(7 -ii:)Eso 


(A.5) 




= 1.06(7 -/:)eso 


(A.6) 


(7 - A')ciT 


= (7 - K)c - 0.105(7 - K)c - 0.001 


(A.7) 




^ Q.2>95{J-K)c 


(A.8) 




= 0.895 X 1.06(7 -7i:)Eso 


(A.9) 




= 0.949(7 -/s:)eso 


(A. 10) 



This is very similar to the colour equation given by Persson 
et al. (11558 ( (hereinafter P98) (7 - K)cyy = 0.954(7 - K\co 
between the CIT and LCO systems. As the slopes of these two 
relations are almost identical and the colour range is very small 
(» 0.4 mag) we conclude that we only need to apply a sim- 
ple offset to the ESO data to bring the dataset onto the LCO 
instrument system. 

A.2. Transformation to the CIT standard system 

On the nights Dec. 27 and 28, 1998 at LCO we observed eight 
different standard stars from the list of P98 with a small range 
in colour (0.24 < (J - K) < 0.37). The stars were observed 
several times during the nights. 

We have performed synthetic aperture photometry on these 
stars using the photcal package within IRAF and we deter- 
mine the zero point shifts Zj and Zk to the LCO standard sys- 
tem (P98) on the basis of the following relations: 

-/lCO = JLCOinst + kjX+Zj (A. 11) 

-'^LCO - Ki^CO inst + kxX + Zk (A. 12) 

X indicates the airmass. We have adopted the canonical air- 
mass terms of kj - 0.10 and kK - 0.08 as quoted by P98 as our 
data did not suggest a significantly different value. 
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Similarly, synthetic aperture photometry was performed on 
isolated stars in the science frames. Growth curves were fit- 
ted for both the science and standard star data. On the basis of 
these data a number of reference stars in each field could be cal- 
ibrated. The remaining time series photometry, which was al- 
ready transformed to the LCO instrument system as described 
above, could then be offset to match these reference stars using 
the transformations from eauations lA. 1 IIA.'O 

The resulting light curves were averaged for each dithered 
observation set consisting of typically 5 or 6 exposures. The av- 
eraged magnitudes were finally transformed to the CIT system 
using the transformation given by P98: 

^CIT = ^LCO (A. 13) 

{J-K)cn = 0.954(7 -/:)lco + 0.015 (A.14) 

The K and (7 - K) photometric measurements are tabulated 
in Tab^and the K and (7 - K) light curves are plotted together 
with the (V - K) colour curve in Figs lll6l 

The (y - K) colour was derived by combining the V data 
from Gieren et al. (i2000a.) with the K-bwd data presented here. 
As the data were obtained at different times it was necessary to 
first smooth the /T-band data and then interpolate in phase to 
derive a (V - /T) colour corresponding to each V observation. 
The smoothing was done by computing the weighted mean of 
5 points in phase and magnitude and interpolating linearly be- 
tween the resulting points. This approach results in fewer ar- 
tifacts than a Fourier fit and thus gives a better representation 
of the data. Due to the low amplitude of the /T-band light curve 
any remaining inaccuracies in the smoothed /T-band light curve 
causes only slight changes in the resulting (V - K) curve. 

The photometric errors as returned by DoPHOT are tab- 
ulated with the magnitudes in Tab^ Due to the fact that the 
stars have been observed with individual pointings as they did 
not fit into the limited field of view of the near-IR cameras, the 
zero points for each field have been determined independently 
using non-variable stars within the fields. We estimate that the 
uncertainty on these zero points is 0.04 mag. 
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